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The explicit formulae for radiative correction (RC) cal- 
culation for elastic ep-scattering is presented. Two typical 
measurements of polarization observables such as beam-target 
asymmetry or recoil proton polarization, are considered. Pos- 
sibilities to take into account realistic experimental accep- 
tances are discussed. The Fortran code MASCARAD for pro- 
viding the RC procedure is presented. Numerical analysis is 
done for kinematical conditions of TJNAF. 



I. INTRODUCTION 

Precise polarization measurements of nucleon form 
factors in electron scattering is an essential compo- 
nent of new-generation electron accelerators such as CE- 
BAF. This unprecedented precision requires knowledge 
of higher-order electromagnetic effects at a per-cent level. 
The purpose of our work is to analyze radiative correc- 
tions in elastic electron proton scattering and develop 
proper computational techniques that could be used in 
experiments at Jefferson and other electron accelerator 
labs. 

The large statistics obtained in the experiments at CE- 
BAF 1^ requires adequate taking into account radiative 
effects. 

Both spin averaged and spin dependent parts of the 
cross section (cr" and a^) can be presented as 



^obs ^ 



u.p 



u,p 



(1) 



Both the factorized correction 6 and unfactorized cross 
section coming from bremsstrahlung process contribute 
to cross section. When polarization asymmetries of elas- 
tic processes are considered, the factorized part of total 
RC has tendency to cancel but unfactorized part can give 
important contribution. This is the reason why we pay 
special attention to the effects which have different con- 
tributions to polarized and unpolarized parts of elastic 
cross section, namely 

• A contribution of hard radiation. Its contribution is 
restricted by using of kinematical cuts on deviation 
from energy-momentum conservation. 

• A higher order correction, which contribute to the 
polarization asymmetry (due to unfactorized ef- 
fects) even at the leading order. 



• A contribution of box-type diagrams which are a 
typical example of the unfactorized correction. 

First and second correction are so-called model in- 
dependent RC. They do not depend on any assump- 
tions on nuclear interactions. In our case it means that 
model independent correction can be written in terms 
of elastic nucleon formfactors. The third correction is 
model dependent because it requires the consideration 
of two-photon exchange graphs, which cannot be calcu- 
lated without any additional assumptions. Its calculation 
in the soft photon approximation, which gives identical 
corrections to both polarized and unpolarized parts of 
the cross sections, was done in refs. ||,^. 

A modern approach of solving the task on RC calcu- 
lation assumes the exact calculation of the lowest order 
model independent correction. The task can be solved 
exactly. Uncertainties of the model independent RC can 
come only from fits and models used for structure func- 
tions. The second order model independent correction 
also can be performed exactly, however usually it leads 
to cumbersome or lengthy formulae, and realistic approx- 
imations are preferred. In this case there is an additional 
uncertainty coming from using approximate formulae in- 
stead of exact ones. The calculation of model dependent 
correction requires additional assumptions about hadron 
interaction, so it has additional pure theoretical uncer- 
tainties, which are hardly controlled. So usually only 
those model dependent correction are considered, which 
are expected to strongly influence extracted observable 
quantities. In our case such a correction is box-type 
graphs. 

There are two basic methods of calculation of model 
independent QED radiative correction. First one is con- 
nected with introducing of an artificial parameter sep- 
arating the momentum phase space into soft and hard 
parts. One can find a classical review introducing this 
formalism in ref . . However the presence of the artifi- 
cial parameter is a disadvantage of this method. For the 
soft photon part the calculation is performed in soft pho- 
ton approximation, when its energy is considered to be 
small with respect all momenta and masses in the task. 
So this parameter should be chosen as small as possible 
to reduce the region evaluated approximately. However 
it cannot be chosen too small because numerical insta- 
bilities of hard part calculation can occur. In the end of 
seventies Bardin and Shumeiko developed an approach 
M of extraction and cancellation of infrared divergence 
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without introducing this artificial parameter Q Later a 
lot of calculations were performed within this approach 
and a lot of FORTRAN codes created to deal with numer- 
ical calculations. The most known of them are TERAD 
and POLRAD. The detailed review of the approach is 
presented in ref. Q. In this paper we use this approach 
to calculate the RC of the lowest order to transferred po- 
larization and asymmetry in elastic electron-proton scat- 
tering. The method allows to calculate the correction 
exactly. However we will neglect the lepton mass m ev- 
erywhere where it is possible. 

The observed cross section of the process 



e{ki) + N{p)~^e'{k2) + N{p2), 



(2) 



is described by one non-trivial variable, which is usually 
chosen to be squared of momentum transferred. There 
are two ways to reconstruct the variables, when both 
lepton and nucleon final momenta are measured. In the 
first case it will be denoted as = — (fci — fe2)^, and for 
second case it is Q\ = —{p2—p)^- It is clear that there is 
no difference between these definition at the Born level. 
However at the level of RC the symmetry is breaking 
because the radiation of the lepton is considered only. 
We consider both cases in this paper. In the first case 
the structure of bremsstrahlung cross section looks like 



da 
dOf 



(3) 



where /C is a kinematical coefficient calculatable exactly 
in the lowest order. It depends on photon variables. JF^ 
is a bilinear combination of nucleon formfactors depen- 
dent on Q\ only, which is function of photon momentum. 
Usually only final momenta are measured in the part of 
the space. It is controlled by function of acceptance A, 
which is 1 or in depending whether the final particles 
make it in detectors or not. The integral ^ should not 
be analytically calculated for two reasons. The first one 
is dependence of formfactors on Q\. We avoid to use 
some specific model for them. The second one is accep- 
tance usually very complicated function of kinematical 
variable dependent on photon momentum. 

For the second version of reconstruction of transfer mo- 
mentum squared the structure of the cross section is 



da 



(4) 



In this case the squared formfactor is not dependent on 
photon momentum and for An kinematics {A — 1) this 



^Toll for that is non-positively defined hard photon contri- 
bution, which can lead to difiiculties for MC generation see 
. In this paper detailed comparison of explicit formulae ob- 
tained within two considered approaches is given for the case 
of deep inelastic scattering. 



integral can be calculated analytically. In the experimen- 
tal conditions at JLab |^ , both final electron and proton 
were detected in order to reduce background. However 
elastic scattering kinematics was restored by the final 
proton kinematics, while electron momentum was inte- 
grated over. Therefore, formalism related with Eq.(U) 
applies for this case. 

The calculation of RC to the asymmetry in measure- 
ment of recoil proton for leading and next-to-leading lev- 
els was done within the method of structure functions [|| . 
This method was chosen because it allows to take into 
account the higher order correction simultaneously with 
enough accuracy of the lowest order RC. However only 
the case of measurement of final lepton was considered 
there. The traditional calculation based on a covariant 
method of infrared divergence cancellation is performed 
for the two experimental measurements of polarization 
observables: 

• measurement of polarization asymmetry within 
leptonic variables; 

• measurement of asymmetry in recoil proton within 
hadronic variables. 

This approach allows to take into account the lowest or- 
der RC exactly and to calculate RC within experimen- 
tal cuts. The treatment of the box-type correction and 
their influence on transferred polarization requires differ- 
ent methods and will be a subject of a separate investi- 
gation. 



II. KINEMATICS AND BORN PROCESS 

The Born cross section of the process can be written 
in the form 



dao 



Apk 



167rS'2 ' 



(5) 



where S = 2fcip. Kinematical limits for are defined 
as 



< < 



A., 



, X,^S^- Am^M^, (6) 



where (m, M are the electron and proton masses). Due 
to axial symmetry the integration over azimuthal angle 
4> can be performed analytically. However in our case 
kinematical cuts are dependent on this angle so we will 
consider two-dimensional Born cross section 



][,f2 ^n2 

dao = T^dTo = Af, 
Apki 



2 dQ' 
° 327r2S'2 



Born matrix element is 



(7) 



(8) 
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We use standard definitions for (unpolarized) leptonic 
tensor and for hadronic tensor 



with 



w 



(9) 



(10) 



and T.p = 

= ATpM^Gl,, T2 = 4M2^l^±^i^ (11) 

It is convenient to define the contractions 

W\ = L%wl, ^ 2Q\ (12) 
2^?l = Ll,,wl^ = j^{S{S ~ _ M^Q^). (13) 

As a result for the Born cross section we have well 
known formula 



(14) 



(15) 



which is reduced to 

d(7Q ^ A'Ka'^ G\ + TpG\j 
dQ2 ~ -Qi i + rp 

in ultrarelativistic approximation. 

A. Polarized part of cross section 

We consider two polarization measurements: 



1. Initial proton is polarized and final electron is mea- 
sured to reconstruct Q^. In this case there are four 
experimental situations for asymmetry definition: 
target is polarized along (perpendicular) to beam 
or q {q = p2 — p). Corresponding polarization 4- 
vectors are denoted as rjL (iit) or 77' (?7y). 

2. Polarization and momentum of final proton are 
measured. Two polarization states should be con- 
sidered: final proton is polarized along (77^) and 
perpendicular (77^) to q. 

If the polarization vector is kept in a general form the 
polarization part of hadronic vector can be presented in 
the form 



(16) 



i=3 



with 



M 



W^,^ = ie,,^Xa ^^3 ■ (17) 



For the consideration of the case of initial polarized par- 
ticles we have to choose the corresponding representation 
for the polarization vector and structure functions in the 
form 



^3 — —2AI^GeGm, 



= -M^Gm 



Ge — G 



M 



(18) 



when the total hadronic tensor — W^^^ + W^^^- In 
the case of consideration of final polarization states (77 — s- 
77') the same formulae for structure functions (|l8|) but 
different sign for the last one {J- a —Ti)- 

The 4-vector polarization 77 can be expressed in terms 
of 4-momentum of particles in reaction. The four con- 
sidered cases correspond to four representations of polar- 
izations vectors 



Vt 



Vt 



1 



(19) 



1 



{-SX + 2M^Q^ + Am^M^))ki + 



2Af (fci - fca) - p 



1 



(2M 



M 
S^X)ki 4 



+ {2M'Q' + S^S)k2 - Q\S + X)p 



where A = SXQ"^ ~m^\- APQ'^ and A, = Sl+AM'^Q'^. 
We note that the task of the calculation is reduced to con- 
traction of leptonic tensors at the born and RC level with 
uijj'^ using correspondent polarization vector representa- 
tion of general form 



77 = 2(a,,/ci + brik2 + c,jp) 



(20) 



and forthcoming integration. Also keeping in mind the 
case of RC we held the variable X, which is S* — at 
the Born level and depends on inelasticity S — — v for 
RC; Sx = S ^ X. The definition of inelasticity v is given 
in the next section. It should be noted that we do not 
consider effects of normal polarization. It allows us to 
keep only three basis vectors in Eq.(pO[). 

Polarization part of the Born cross section is given by 
Eq.(p^ for two additional terms in the sum over 7=3,4. 
Functions 6^ have the form 



qB 
^3 



^4 



2777, 

mQ^qrj 



IvP 



{qr, k2^-ri^Q^), 

m-k20, 



(21) 
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where lepton polarization vector can be taken as 
2 



m p) 



(22) 



If the is calculated in terms of hadronic variables 
the polarization vector expansion looks like this 



'nL,T = '^ia'L,Th + bL,T{P - P2) + c'l^tP) 



(23) 



with 



0, b'^ 



Q2 + 



M 



2MJ\ 



M 



and 



Q^S + 2APQI 

2\/A/i\/Am 
QH2S~Ql) 



(25) 
(26) 
(27) 

M 

Here A,, = SQ^{S - Ql) - M^Qi - tu'^Xm and Am = 
Q^{Q^ + 4M^). Similarly to the case above we hold the 
variable Q^, which is and + u for Born and RC 
cases respectively. The quantity u is related with invari- 
ant mass of unobserved state. It is also called inelasticity 
and is defined below (for the case of hadronic variables 
see (H)). 

It easy to verify that 4-vector 77^ and 77^ satisfy the 
necessary conditions: 



2Va,,va 



Vt 



-1. 



In the rest frame of the scattered electron p2 = {M, 0) 
the vector of longitudinal polarization looks as 



(0,n), n^ = l 



(28) 



The direction of 3-vector n coincides with the direction 
of 3-vector p2 in lab. system. Therefore, 77^ indeed de- 
scribes the longitudinal polarization of the scattered pro- 
ton. 4-vector rjip has the form (0,TO),n • m = in both 
lab. system and rest frame system of the scattered pro- 
ton. Thus, it describes the transverse polarization in the 
scattered plane. It can be defined up to sign only. 

In the case of longitudinal polarization contributes 
the term accompanied with multiplier G\.j in the spin- 
dependent part of the cross-section, because the part, 



that is proportional to GeGm goes to zero for rj = rj 



II. 

The situation is just contrary in the case of transverse 
polarization. The simple calculation gives: 



Gm 
~G~ 



kip + kiP2 



E V M2 ^4{kip + kiP2Y + q^{AAP - q2) 



(29) 



It easy to verify that in the Breit system, where p = 
{E,-^2), p2 = {E,q/2), q = (0,g), the right side of 
Eq.(p9|) coincides (up to sign) with the expression given 

in Q. Indeed, in this system ei = £2 = 2singfl/2 (^b is 
the electron scattering angle in the Breit system) and 



ikip + kip2) 



1 



y^4{kip + kiP2Y + q^iU'P - g2) 2cos6ls/2 ' 

III. RADIATIVE EFFECTS 

A. Leptonic variables 

For cross section of radiative process 

e(fci) + N{p) e'{k2) + 7(fc) + ^fe), (30) 
we have again 



M2 

ddr = - — ; — dVr- 
Apki 



(31) 



The cross section of the process depends on Qf which for 
simplicity is referred in this section as Q^. Phase space 



dTr 



1 d^p2 d^k2 d^k 
(27r)5 2p20 2fc20 2fco 



5{p + ki 



k~P2) (32) 



can be parameterized in terms of three variables: inelas- 
ticity V = A? — AP (A = p + fci — ^2), T — kq/kp and 
angle 1/)^ between planes (q, k) and (ki, Using the 
result (V.7.7) of we have 



dVr 



dQ^ 



A{2nYS J 4VA 



dv 



dr- 



V 



l2^2 , A ]\/r2n2\ 



d^k, (33) 



where (A, = (v + Q^Y + ^A/PQ^) we use a variable r in- 
stead of standard t 



t = Q^ + v~R- 



R = 



l+T 



(34) 



It allows to present final result in the form close to circle 
of papers |0-fi3. Limits are defined as 



1 



(35) 
(36) 

(37) 
and 



2m2 

2Q^Xs - QHS + + AP)) 
Q^S + 2m2Q2 + VAIVA;^ 



2MV, 



max.min 



±a/a: 



(38) 



Matrix element squared of radiated process is 



4 



(39) 



The leptonic tensor of the radiative process is standard 
and can be found, for example, in We note that 

here we use more standard definition of the tensors (with 
additional factor 2 comparing to ones from the paper). 
For contractions we have 



(40) 



Functions 6ij are similar to ones given in Appendix B of 
rcf. ||T^. However there they are integrated over 0^.. We 
keep this integration because of possible dependence of 
acceptance function on this angle. Explicit form of the 
functions in our general case is discussed in Appendix. 

We note that the well-known formula of soft photon 
approximation is immediately obtained when keeping the 
term with j = 1 and restricting integration over (small 
photon energy) v as vi < v < V2 <^ S (small photon 
energy): 



dar 



2a 



dQ2 ~ TT ^^"^ 



l)log 



V2 dap 
'^dQ^ 



(41) 



where l„i = \og{Q^Jm^) . For angular integration the 
formula (27) of ref [IJ] was used. 

Straightforward integration over photon phase space 
is not possible because of infrared divergence. The first 
step of the treatment is the identity transformation of 
the integrand 



(42) 



where ap finite for /c — > (here and below we use 
the short notation for differential cross sections: an = 
daii/dQ^ and so on). There is some arbitrariness in def- 
inition of aiR. Only asymptotical expression in the limit 
k ^ is fixed^. In our case we construct aju using 
the term with j = 1 in (40) and formfactors estimated 
in Born point. This term is factorizcd in front of Born 
cross section as 



(Tq- 



2ka 



IR, 



IR 



ki k2 
2kik ~ 2k2k 



(43) 



As a result infrared part can be written in the factorizcd 
form 



TT — 



a 
-( 

TT 



(44) 



There is once more limitation. We must provide the condi- 
tions of applicability of the theorem about changing of order 
integration and limitation. For example the uniform conver- 
gence is required. Practically it means that we may subtract 
the quantity with the same denominator 



The quantities 63 and Sh appear after additional split- 
ting the integration region over inelasticity v by the in- 
finitesimal parameter v 



Ss = — / dv 



Jn-l I 



(27r/i)"-4A:o 



Sh = — / dv 



f^FinSHA-kr-M'). (45) 
ko 



We note that in contrast to Mo and Tsai formalism this 
artificial parameter v completely cancels in final expres- 
sions. The way to calculate these integral was offered in 
Q (see also [|4j and review Q). In our case we have 



2 P, 



V 



(46) 



These contributions have to be added with result for ver- 
tex function which is standard 



-2{ PiR 

1 



log — ) {Ir, 



1) 



(47) 



2 2 " 



2 + — 



For this sum we have 

a 

TT 

where 



Ss + Sh + Sv 



^VR, 



Sinf 


--(I - 


7T 


SvR 








~ 




+U2 (1 



1 log 



(48) 



(49) 



S 



S{S~Q^)) ~ "6 

Here we used ultrarelativistic expression for function 
from 

Finally the cross section that takes into account radia- 
tive effects can be written as 

a-obs = CToe''"'^ (1 + SvR + Syac) + crp- (50) 

Here the corrections Sinf and Syac come from radiation 
of soft photons and effects of vacuum polarization. The 
correction Svr. is an infrared-free sum of factorizcd parts 
of real and virtual photon radiation, ap is the infrared- 
free contribution of bremsstrahlung process: 

f , f dr 

252 



dv 



1 



9fe X 



(51) 



E 



4F 



IR^i 



. ^j = l -^h RQi 

Here A is integrated over the cj) acceptance function 
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B. Hadronic variables 

For the cross section of radiative process 
e(fci) + Nip) e'{k2) + i{k) + N{p2) 
we have again 

dcjr — — ; — dVj. 



(52) 



(53) 



The parameterization of photonic phase space, inte- 
gration over it developed for so-called hadronic emission 
within considered Bardin and Shumeiko approach can 
be applied to this case, if we picture the corresponding 
Feynman graph upside down. Thus phase space (|3^) can 
be parameterized in terms of three invariant variables 
0,0: inelasticity u = {k2+k)^-m'^ = 2fc2fc, w = 2kik 
and z — 2p2k 



dTr — 



dQ^d<t) 
(47r)4S' 



du 



dw 



dz 



(54) 



where Rz comes from Gram determinant {16Rz = 
A{ki,p,p2,k)) and coincides with standard /J^-function 
appearing in Bardin-Shumeiko approach dJl^- Explic- 
itly 

Rz^ Azz^ ^2BzZ + Cz. (55) 
For completeness we give coefficients in our notatiori|^ 



Bz = u{u + Q2)s, - (u - Q'^)Sw - 2m^Q^{u - w), 
Cz = {usq - Sw^ + AM'^uw{Q'^ +u~w) 



(56) 



where Sn 



S - Q'^ + w - u. 



We note that we introduced invariant variable z which 
corresponds to azimuthal angle in (js^) . It is more conve- 
nient for introducing explicit expressions for experimen- 
tal cuts. 

Limits are defined as 



1 



2M2 

2Q^{Xs-Q^{S + m^ + AP)) 
Q^S + 2Af2Q2 + 



M 



and 



2{u + m?) 



{Q"^ +U + 2m^ ± ^/Xq 



(57) 
(58) 

(59) 



^Q^ in this subsection is Q\ 



Limits Zrain max are defined as solutions of the equation 
i?. = 0. 



Bz±Vd, 



(60) 



where 



D = iiM^Xg + m^Q^ - Q^S{S - ~ u)) x 

x(m^(w - u)^ + uw(w - u - Q^)). (61) 

Two solutions of the equation Z) = give limits on w 

(!)• 

Matrix element squared of radiated process is calcu- 
lated as 



Ml 



where 



Q 



jLl,W^. = -^{TiR + T; + n + ni) (62) 



T. = 



MS* 



2^(Q z - 2Q^Su + Q^uz - u^S) 



u w 



2Q'a'L^T^i (2 



-{Q'^Su - Q^uz - 2S^u + Su^ 



S . 



T34 = 



+2Suz - 2uz^) + -i-2Q^ + AQ^S - 2Q^u 
w 

-2Q^z + 2Su - u^) - 2S^ + Sw + 2Sz^ , 



X 2 — {Q'^z -Su~ 2Sz + 2z^ 



25zQ2 



--(25 -w^2z) + -(2Q2 - 25 + m) 
u w 



The contribution Tm is 

To = -^{2M^Wl^tQ^ - ^l^tQ^ + 2c'l^tS) 

+Q2(q2 _ 2S)Ti{a'L^T ~ ^b'^.T - c'l,t)) ■ (63) 
Structure functions J-'j and J-a are defined in (nSf) and 



are functions of . We note that Born cross section can 
be written in terms of Tq for a^,, 6^, c,, taken for u ^ 



(To = 



1 n 



The radiative cross section can be given as 
—a^ f dudwdz 



(64) 



0-7? = 



-.A{TiR + n + Ti + n^). (65) 
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Performing the explicit integration over the region of 
smah energies (0 < ui < U2 ^ m'^ , M'^ , Q"^ , S) we ob- 
tain the weU-known result in the soft photon limit: 



du 



dw 



dz 



(jR = (Jo / du 

TT 



1 



u u{u + Tin?) 



Q 1 '^max 

-^=log 



(72) 
(73) 



2a U2 , , 

log [Im - IjO-Q- 

TT Ui 



(66) 



The radiative cross section has infrared divergence, so 
for this case the identity transformation like (E3) has to 
be performed also. However formfactor is not depen- 
dent on photon variables, so only the acceptance func- 
tion should be subtracted. Corresponding integrals looks 
similar. As a result infrared part can be written in the 
factorized form 



4Q4^2 J J ~~ J 

Wmin Zmin 

A{Tjn + T3 + T4 + T34) - tPr 
with TPj, = TPj,{To ^ To^) 

C. Kinematical cuts 



In this section we show how experimental cuts can be 
introduced to this approach. As an example we consider 
conditions of experiment at CEBAF. 



oiR = —Oj.' cTo = -((Js + Oh)<7o- (67) 

TT TT 


RECOIL PROTON ^ 






The quantities Ss and Sh appear after additional split- 
ting the integration region over inelasticity v by the in- 
finitesimal parameter v 














BEAM DIRECTION 



6h = 



du 



ko 



(68) 



where Ah = fci + P — P2 and Fjj^ is defined in eq.(43) 
The integration gives the following explicit results 



6^ = 2( P^^ 



5'k = 2{l^ - 1) \o; 



log (?,„ -!) + ! + U - ll, - — , 

m/i / D 

1 



2 

^2 



log- 



^Iwilm + ^ii> ^ Iv) Li2 I —r^ ] , 



1 1 """^ 



(69) 



Sum of them and the contribution of vertex function gives 
again infrared free result 



-Li2 



(70) 



e 



FIG. 1. Recoil proton angles definitions. The beam and 
the hadron-arm detector (HRS) define the horizontal plane. 

The following restriction on momentum of final elec- 
tron and proton have to be done. We consider high- 
resolution spectrometer (HRS) |l| as a rectangular area 
with some energy acceptance. We describe this rectan- 
gle by two angles between corresponding planes: 0^''' and 
9y'P. For one of the scatteredparticle these angular def- 
initions are given in Figure The upper index corre- 
sponds to detected particle: electron or proton. The 
angles and momenta of final particles being measured 
in lab. frame have to be expressed in terms of kinemat- 
ical invariants. The simplest way to do it is to apply 
a formalism of Gramm determinants, which detailed de- 
scription can be found in Ref . (l^ . We will give starting 
expressions in terms of four-momenta and Gramm deter- 
minants as well as explicit results for the used (in (jsj)) 
invariant variables. In terms of Gramm determinants the 
momenta (in the rest frame of p) of final proton and elec- 
tron are given by the formula 



The cross section that takes into account radiative ef- 
fects can be written as 



IP2I 



MP,P2) |, ,2 A(p, fc2) 
2 |k;2| = 5 • 



C^obs — CToll + Syjf + Svac) + CTp. 

The explicit expression for ap is 



It gives immediately for the Born process 



(71) 



IP2I 



2M ' 



Iksl- 



2M 



(74) 



(75) 
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both for Born and radiative process. Here Xz = {S — 
Ql - zif - iM^m\ 

Cosine of polar angle of the direction of final proton 
(with respect to the beam direction) is defined as 



C0S^p2 = 



Aip, fci)A(p,p2) 



It gives 



sm 



7p2 



(76) 



(77) 



where A^^^ SXQ"^ - S^M^. We have to use X = S*- 
and 5^; = for Born process X = S — — u and 
Sx = + u for radiative one. In terms of the angle the 
horizontal and vertical angles of proton momentum are 



sm sm ( 



p2, 



tan 0- 



p2, 



(78) 



where </) was introduced in (0). 

At the Born level all vectors (p2, ki and k2) are in the 
same plane. However for the unobserved photon there is 
a nonzero angle between planes (ki, P2) and (ki, k2). 
In terms of Gramm determinants it is defined as 



cos A(/) : 



Q f P ki P2 
p ki k2 



A{p, ki,k2)A{p, ki,p2) 



(79) 



Explicitly we have 



sin A(j) ~ — 



■ [{S{u — w) — ziu + Q^wY 



+AQlM^uw + Q^{2w + zi)zi(Q2 + „) 
+zluQ'^ + 2ziQ'^uw -2Szi{u + w)Q^] (80) 

where A,, = SSq, - 2KPQl, Sg, ^ S - - zi, Ql = 
+ U — W and zi = z + u — w. Now we can define angles 
of final electron for radiative process 

sin 61^ = sin(0 + Ac/)) sm6k2, 
tan 9"^ = cos((/) + A0) tan 6*^2, (81) 

where 6^.2 is polar angle of final electron. For it we have 



.2^_ _4g2A,,M2 



sm = 



q2 02 



(82) 



the observable quantities in polarized ep scattering mea- 
surements. The numerical results of applying this code 
for k inema tica l con ditions of Jlab are given in Subsec- 



tions [VB 



and [V C with leptonic and hadronic variables 
respectively. In last Subsection we discuss the influence 
of experimental cuts on observables quantities in polar- 
ized scattering. 



A. FORTRAN code MASCARAD 

There are two variants of the code: mascaradJ.f and 
mascarad_h.f dealing with leptonic and hadronic vari- 
ables respectively. First code does not require any exter- 
nal libraries. However the histogramming by HBOOK 
can be optionally included in second variant. In this case 
mascarad_h.f requires CERNLIB installed. In exter- 
nal file one can choose kinematical variables, accuracy of 
calculation and value of a cut on inelasticity. An option 
to include kinematical cuts described in Section 
available for mascarad_h.f. As an output one 



IHC 



las the 

value of the Born cross section, radiative correction factor 
(with estimation of statistical error) in chosen kinemati- 
cal points. 



B. Numerical results. Leptonic variables 

Absolute and relative corrections to asymmetry can be 
defined as (see (0)) 



AA, 
A,: 



A, 
A, 



A;0 

-Am 



(l + J)ag + a^ 



A,o 



1 + S + Su 



(83) 



(84) 



where index i runs over all considered cases: i ~ 
L,T,qL,qT] Su,p — ^i^'^/cg'^. Here the correction S 
is usually large because of contributions of leading log- 
arithms. However it exactly canceled in numerator of 
expression for correction to asymmetry. That is a rea- 
son why the correction to cross section can be large, and 
correction to asymmetry is relatively small. 



IV. NUMERICAL RESULTS 



In this section we presen t the FORTRAN code MAS- 
CARAD (Subsection IV A) developed on the basis of for- 
malism presented in the last sections. This code uses 
Monte-Carlo methods to calculate radiative correction to 
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FIG. 2. Born (solid line) and observed (dashed line) asym- 
metries vs . No kinematical cuts on inelasticity ware used. 
E = A GeV. 

Born and observed asymmetries are presented in Fig- 
ure ^ Four lines corresponds to four considered cases de- 
fined in sections If A. No cuts were used for the missing 
mass. As a result hard photon contributions gives dif- 
ferent contributions to spin averaged and spin dependent 
part of cross sections due to unfactorizing properties. For 
longitudinal asymmetries 5p > Su and there are positive 
contributions to RC. Opposite situation with transverse 
asymmetries. 

One can see that the transverse asymmetry At in re- 
spect to beam direction has large correction. It is not in 
contradiction with other plots in this figure. The polar- 
ized part of cross sections in the cases (L,T) and (qL, qT) 
are related to each other by some unitary transformation 



_dal_ 

dQ^dv 



cos 7 



dQ^dv 



sm7 



dQ^dv ■ 



(85) 



Because of dependence of polarization vectors on inelas- 
ticity (see ( |l9| ) and definitions there) this angle 7 is a 
function of v 



sin^ 7 



4M2A 



COS 7 



SS^ + 2M2g2 



•\/ As Aq 



(86) 



and only unintegrated cross sections are related as (85). 
This sine hardly suppresses the cross section of the hard 
photon. Weighed with the sine and cosine, cross sections 
in r.h.s of (^5|) have the same signs, similar magnitude 
and compensate each other. As a result 6p <C <5„ and the 
asymmetry At has a large negative contribution. 

In practice RC to the asymmetries can be essentially 
reduced by applying the cut on missing mass or inelas- 
ticity which is also measured quantity in elastic electron- 
proton scattering. In Figure p it is shown how these rel- 



ative corrections depend on value of the cut on missing 
mass or inelasticity. 




FIG. 3. Relative RC to asymmetries defined in ( |83| ) vs 
value of inelasticity cut for Q^=3 GeV^; E = 4 GeV. 



C. Numerical results. Hadronic variables 

Similarly to the case of leptonic variables let us define 
relative RC to the ratio of recoil nucleon polarizations 
Pt/Pl as 



A = 



Pt/Pl -P^/ Pi 
P^/Pl 



(87) 



In Figure ^ this correction is given for several values of 
cut on missing mass. It can be seen that radiative cor- 
rection in this case is smaller then in the case of leptonic 
variables. 

However in practice experimental situation can be 
complicated then applying single kinematical cut on 
missing mass. In the last and future experiments ppg|] 
at JLab on measurement ratio of elastic formfactors of 
proton the all events appearing in detectors^ were ac- 
cepted for analysis. To calculate RC in this situation we 
have to apply all cuts discussed in previous section. The 
Table | gives the results for the last |0] and future |T9| ] 
experiments (upper and below the intersection line). 



*Here we do not consider possible losing events in detector. 
By the other words we take into account only geometrical 
acceptance. However influence of apparatus acceptance also 
can be taken into account in our approach by introducing 
some map of device acceptance to MASCARAD. There is a 
certain place in the code to do it. 
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As we can see the RC does not exceed 1%. 



V. DISCUSSION AND CONCLUSION 





-0.002 
-0.004 
-0.006 
-0.008 
-0.01 
-0.012 
-0.014 
-0.016 
-0.018 



1 - v_m=0.1 GeV^ 

2 - v_m=0.3 GeV^ 

3 - no cut on V m 



5 

Q^ GeV^ 



FIG. 4. Relative RC to the ratio of recoil proton polariza- 
tions vs for three values of inelasticity cut. 



TABLE I. 


The results for RC to 


asymmetry Pt / Pl 


E,GeV 


Q^GeV^ 


A % 


0.934 


0.45-0.53 


-1.01 


0.934 


0.77-0.81 


-1.54 


1.821 


1.11-1.25 


-0.95 


3.395 


1.37-1.59 


-0.59 


3.395 


1.65-1.89 


-0.64 


4.087 


1.75-2.01 


-0.62 


4.090 


2.30-2.64 


-0.71 


4.087 


2.77-3.17 


-0.80 


4.090 


3.27-3.67 


-0.95 


4.845 


3.5 


-0.80 


4.845 


4.2 


-0.96 


5.545 


4.9 


-0.95 


6.045 


5.6 


-0.97 



In this paper we applied the approach of Bardin and 
Shumeiko for calculation of the model independent radia- 
tive correction of the lowest order in processes of elastic 
electron-proton scattering. Current experiments on the 
process measure different polarization observables such 
as spin asymmetries in different configurations of polar- 
izations of initial particles and the ratio of recoil proton 
polarizations allowing to access the ratio of electromag- 
netic formfactors of the proton. That is way special at- 
tention was paid to radiative correction to polarized parts 
of the cross section. 

The chosen method of calculation allowed to obtain 
explicit formulae in the cases of so-called electron and 
hadron variables. They correspond to the cases when 
kinematics of the measured process is reconstructed from 
the momentum of final electron and proton respectively. 
It was shown that in spite of the formulae for Born case 
are exactly the same for both cases, all ingredients of the 
RC calculation are different in these cases. Physical (or 
kinematical) reason of it is the fact that in the first case 
the radiating particle (electron) is measured, but in the 
second case it escapes unmeasured or is integrated over. 

Explicit formulae for the lowest order radiative correc- 
tion are exact if not to account ultrarelativistic approxi- 
mation taking off the terms ~ vr? jQ^ . However it is not 
a limitation of the approach. These terms can be easily 
restored in the case of need (for example, if they will ap- 



ply to to muon scattering and the accuracy of ■ 



will be required. 

Contrary to the case of inclusive DIS the integration 
here is left for numerical analysis. All integrals are finite 
after used the procedure of covariant cancellation of in- 
frared divergences. This form of answer allows to include 
acceptance effects to the integrand. The function of ac- 
ceptance usually depends on final angles and momenta, 
which can be expressed in terms of integration vari ables . 
Corresponding way to do it is discussed in Section III C . 

On the basis of the exact formulae the FORTRAN 
package MASCARAD was developed. It includes codes 
both for the electron and hadron variable measurement. 
The applying of the package to the radiative correction 
procedure allows to provide the account of the model in- 
dependent correction in current measurements (including 
polarization) of elastic electron-nucleon scattering. Our 
numerical analysis shows that radiative effects can be im- 
portant especially in the cases of transversely polarized 
targets. However using kinematical cuts such as single 
cut on inelasticity or a cut on kinematical variables of 
second (undetected) particle allows to reduce the effect 
essentially. 
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APPENDIX 

The functions 9ij are defined the same way as in Ap- 
pendix B of ref. or with more details in in Appendix 
B of ref. [|2j . However there the formulae are integrated 
over photon azimuthal angle (or equivalently over z, see 
(B7) of ref. Below we define the procedure how the 

explicit form of this functions can be written in out case. 

All formulae (B1-B2) of ref. flltl or (B.l-B.ll) of ref. 
1^ can be kept for our case. Instead of functions F's 
from (B5) or (B.12) ^ we use the following expres- 
sions 



F 



F F 

, Fi+ = 1 ,F2±=F{ ± 



Zl 



zi 
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